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Abstract 

We study the vacuum stability of the singlet Majoron model using full renor- 
malization group improved scalar potential and Monte Carlo techniques. We show 
that in the perturbative regime of the various free parameters, the vacuum stability 
requirement together with LEP limits is passed by 18% of the parameter space if 
the scale of new physics is 10 TeV and 6% if the scale is 10 14 GeV. Moreover, if the 
baryogenesis condition for scalar couplings is required, no portion of the parameter 
space survives. 



The singlet Majoron model as one of the simplest extensions of the Standard 
Model (SM) offers an explanation to many problems that remain open in the mini- 
mal Standard Model. In addition to the usual doublet scalar field of the Standard 
Model the singlet Majoron model contains a complex, electroweak singlet scalar 
field. Also, right-handed electroweak singlet neutrinos are introduced to the model 
so, that an extra, global U(l) symmetry appears. This global symmetry is bro- 
ken approximately at the electroweak scale manifesting that the other of the new 
scalar degrees of freedom becomes massive whereas the other remains as a massless 
Goldstone boson. The non-zero vacuum expectation value of the singlet field then 
implies a mass term for the right-handed neutrinos. 

A strong pro in favour to such an extended model of the Standard Model is 
that the weak scale baryogenesis, which appears to be somewhat problematic in the 
Standard Model [Q, |3|, is easier to realize in the singlet Majoron model 0, 0. A 
problem with the Standard Model is that the experimental lower bound for the Higgs 
boson mass ||, ran > 60 GeV, is too high compared to its theoretical upper bound in 
order to avoid the erasure of the baryon number by sphaleron mediated transitions 
P|. In the Majoron model these problems are circumvented by introducing the new 
scalar degrees of freedom ||[ ||. Moreover, the singlet Majoron model together with 
the sew-saw mechanism offers an explanation to the vanishingly small (left-handed) 
neutrino masses. By-products of the sew-saw mechanism, the heavy right-handed 
neutrinos may also serve as a source for cosmological baryon asymmetry due to 
their lepton number violating decays which can be converted to baryon number by 
sphalerons ||. For a review of some other baryon number production mechanisms 
applicable in the singlet Majoron model, see JTTi|j . 

To be a realistic model of electroweak scale physics, any model has to have the 
vacuum of the observed universe stable enough. If the stability is not achieved, 
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the theory can not be fully correct. This stability is, however needed only up to 
some scale, like super symmetry scale or GUT scale, where new interactions becomes 
important: new phenomena appears and saves the stability. In this letter we study 
the stability properties of the singlet Majoron model at zero temperature. We work 
using one-loop perturbation theory and make a full renormalization group (RG) 
improved stability analysis of the scalar potential. Such an analysis for the Standard 
Model has first been performed by Flores and Sher and recently by Sher [12 



for new. They found that the stability of the potential requires that the Higgs mass 
rriH > 75 GeV+1.64(m top — 140 GeV). For top quark mass 174 GeV, which we use as 
an example, m# > 131 GeV. Renormalization group improved stability analysis for 
two doublet model has also been carried out fpL3 l . The singlet Majoron model case 



differs from both above mentioned cases, not only because of different scalar content, 
but also because of the inclusion of the right-handed neutrinos. Their unstabilizing 
effect to the potential may be remarkable like the effect of the top quark has in the 
minimal Standard Model. 

The tree-level potential of the singlet Majoron model reads 

V {H, S) = m 2 H \H\ 2 + m 2 s \S\ 2 + 7\H\ 2 \S\ 2 + (3\S\ 4 + \\H\\ (1) 

where, in a spontaneously broken theory, the mass-like parameters m 2 H and m 2 s 
are negative. The couplings (3 and A are positive with r y 2 < AX/3 which guarantees 
that the tree-level potential has a stable, non-trivial minimum. (See Ref. [|J for a 
detailed study of the potential.) The right-handed neutrinos couple to the singlet 
field with Yukawa couplings gi {i = 1, 2, 3) corresponding to the three generations of 
leptons. In practise we, however, consider only the heaviest right-handed neutrino 
and disregard the other two. The moduli of vacuum expectation values \[2f and 
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\/2/ of the scalar fields H and S, respectively, are given by 

-27m| + A/3m 2 H 



P 



7 2 - 4A/3 : 
-27m| + A\m 2 H 



(2) 
(3) 



7 2 - 4A/3 ' 
and the mass eigenstates are 

ml = \f 2 + (3f 2 ± y/D, (4) 

where D = (A/ 2 - /?/ 2 ) 2 + 7 2 / 2 / 2 - 

The full set of the renormalization group equations needed to stability analysis 
calculated using one-loop perturbation theory reads 
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4A 7H + ^[5 + 12A 2 + i 7 2 ], 

4/37s + ^[5' + 10/3 2 + 7 2 ], 

27 [7h + 75] + ^[7 2 + 2 7 /3 + gTA], 



(5) 
(6) 
(7) 
(8) 
(9) 
(10) 

(11) 
(12) 



where a a = <7^/ (47r) for the gauge as well as for the Yukawa couplings. The param- 
eters B and B' are functions of couplings defined by 

~3a 2 + 2aa' + (a') 2 



B 



B' = 



3n z 



-Aa 2 



-Qn 2 



222 
a± + a 2 + a 3 



The anomalous dimensions of the scalars 
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1H = 
IS = 



[—3a — a' + Aa t ], 
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(13) 
(14) 

(15) 
(16) 



are also needed. The dots stand for the derivative with respect to the logarithmic 
scale variable t = ln(/x/M^), where \i is the effective scale. We renormalize the 
couplings as a,(0) = 0.114, a(0) = 0.0335, a'(0) = 0.0102, a t {fi = 2m top ) = 
m to P /(^ n f 2 ) = 0.0790 for m top = 174 GeV. The right-handed neutrino couplings 



are renormalized using [14 



1 Tfl^ 

ai(ji = 2m Ni ) = 7^yt> ( 17 ) 



where M Ui are the physical heavy neutrino masses. The vacuum stability requires 
now |L5| that inequalities f5{t) > 0, X(t) > and 7(t) 2 < Aj3(t)\(t) holds for all t 



up to the given scale of new physics. Note, that the quadratic terms of the scalar 
potential are not important because we are only interested in large scales. 

The large number of degrees of free parameters makes the integration of Eqs. 
(5) - (12) complicated and difficult to visualize in practise. Therefore, to study the 
RG-equations, we perform a Monte Carlo analysis where the initial values of the 
parameters A, (3, 7, / and the largest neutrino Yukawa coupling gy are generated 
randomly at // = Mz- The logarithms of these parameters are taken to be uniformly 
distributed, so that we are able to cover the whole parameter space rather densely. 
As mentioned, we neglect the other right-handed neutrino Yukawa couplings but the 
largest one. This does not invalidate the analysis due to the structure the Yukawa 
couplings emerge in the RG-equations. The ranges of the parameters are 

RT 3 < /3,7< 10"\ 5 x 10" 3 < A < 10" 1 . (18) 

These bounds are due to the requirements that the baryogenesis analysis of || and 
the perturbation theory are applicable. Moreover we choose 

1 GeV < / < 10 TeV, 10~ 3 < g Y < 1, (19) 

where the lower bound for / comes from the experimental upper bound and the 
see-saw mechanism estimate for the mass of the electron neutrino m Ue ~ m 2 e /(gYf)- 



On the other hand / cannot be much larger than / = 246 GeV because the singlet 
scalar would effectively decouple [|J. (This kind of model is, of course, technically 
allowed, but is physically not very meaningful containing a new, ad hoc symmetry 
breaking scale.) The range of gy was taken to be on the perturbative domain but 
not too small for the second order corrections in /3, A or 7 to be effective. 

From the randomly chosen initial values only those are accepted for which the 
two conditions, 4/3A > 7 2 at fi = Mz and there is no Landau singularity below the 
maximum scale studied, holds. Approximately 73% of the generated values pass 
this test when maximum scale is the supersymmetry scale \x = 10 TeV and 71% 
when the maximum scale is the unification scale \x = 10 14 GeV. Furthermore three 
different cuts are applied to limit the parameter space: 

1. It is required that vacuum is (absolutely) stable up to the given maximum 
scale. In the case of supersymmetry scale 24 % of the generated points in the 
parameter space pass this cut and in the case of unification scale 8 % pass the 
cut. 

2. It is required that the scalar masses calculated from the generated point are 
compatible with the model independent LEP lower limit for a scalar particle 
mass. This means (see |3| for details) that 

cos 2 z T(m + )BR + + sin 2 z Y(m_)BR_ < T(60 GeV), (20) 

where T(m) is the decay rate of Z to a fermion pair and a scalar with mass 
in, BR± are the branching ratios of scalars to ordinary fermions and z is the 
angle of rotation from current to mass eigenstates. 30 % of the points pass 
this cut. 

3. Supplementary to the second limit there is another limit from LEP data. The 
latter relies on a event topology analysis of the reaction Z — > Z*H — > l + l~(uv) 
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when the former relies on observing reduced number of the standard fermionic 
decays of this reaction |§. The resulting limit reads 

cos 2 zT(m + ) + sin 2 ^r(m_) < 4 x T(60 GeV). (21) 

Also this cut is passed by 30 % of the points. 

In Figure 1 the resulting projection of the parameter space in the m + -m^ - 
plane is shown after applying three combinations of the above cuts. Note, that 
these masses are not really the physical ones, because their scale dependence is 
not taken into account. As the scalar masses are not too far from the weak scale 
Mz the correction due to the running of the masses may, however, be neglected. 
The parameter space is reduced to 18 % after applying all cuts in the case of the 
supersymmetry scale and to 6 % in the case of the unification scale. The vacuum 
stability bounds A to be larger than 0.034 (0.067) for /i = 10 TeV (10 14 GeV) 
except some very few points which cover less than 0.03% of the parameter space. 
Therefore it can be concluded that vacuum stability requirement and baryogenesis 
bound A < 0.018 [|J are somewhat controversial. Applying the vacuum stability cut 
and the baryogenesis bound for A together leaves zero parameter space. It is also 
noteworth that Monte Carlo analysis shows that most of the values of the neutrino 
masses passed by the vacuum stability and LEP limits are between 20 GeV and some 
hundred GeV's, i.e. not very small nor very large neutrino masses are favoured. 

Our results infer that the vacuum stability and baryogenesis bounds in the sin- 
glet Majoron model can not be fitted together. However, both limits have some 
uncertainty. Because the baryogenesis bound may be affected by non-perturbative 
effects similar to ones in the Standard Model, the upper bound for the doublet 
self-coupling A may be relaxed. Yet, if our vacuum is not absolutely stabile but its 
life-time exceeds the age of the present universe ~ 10 10 yr, the vacuum stability 
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bound will be somewhat relaxed, too. (For a general discussion, see Ref. |H| and 
references therein.) In the Standard Model the possibility of an un-absolute vacuum 
allows, however, only about 10 GeV smaller higgs mass than the absolute stability 



requirement [12"] . It is not probable that this effect would be large either in the sin- 
glet Majoron model. Taking into account the uncertainties of the bounds it might be 
possible that a very thin right-angular zone beginning at (m_,m + ) « (0,90) GeV 
and ending at (m_,m + ) ~ (90,150) GeV survives both baryogenesis and vacuum 
stability bounds, i.e. the gap between A < 0.018 and A > 0.034 can be filled. 
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FIGURE CAPTION 

Figure 1. Distribution of randomly generated points in the m_ — m + -plane after 
applying combinations of cuts, a) Points (26 % ) which fulfill LEP limits (^) and 
(|2T|). b) Points which pass all cuts when the scale of the new physics is \i = 10 TeV. 
c) Same as b) for /i = 10 14 GeV. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9404268vl 
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> 

O 



80 
60 
40 
10 



> • 



• v../ A.V, ... • 




220 

200 
1 80 
1 60 
1 40 



. . . •* • • * . . 




220 
200 
1 80 
1 60 
1 40 
1 20 
1 00 
80 
60 
40 



o 



40 



60 



80 



1 00 



20 



m_ [GeV] 



